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Abstract 

The fluid-structure interaction is studied for a system composed of two coaxial pipes 
in an annular geometry, for both homogeneous isotropic metal pipes and hber-reinforced 
(anisotropic) pipes. Multiple waves, traveling at different speeds and amplitudes, result 
when a projectile impacts on the water filling the annular space between the pipes. In the 
case of carbon fiber-reinforced plastic thin pipes we compute the wavespeeds, the fluid 
pressure and mechanical strains as functions of the fiber winding angle. This generalizes 
the single-pipe analysis of J. H. You, and K. Inaba, Fluid-structure interaction in water- 
filled pipes of anisotropie composite materials, J. FI. Str. 36 (2013). Comparison with 
a set of experimental measurements seems to validate our models and predictions. 

Keywords: Fluid-structure interaction; Water-hammer; homogeneous isotropic piping ma¬ 
terials; carbon-hber reinforced thin plastic tubes. 


1 Introduction 

This article is part of a series of papers m, 0 , 01 , 0 , m devoted to the investigation of 
water-hammer problems in fluid-filled pipes, both from the experimental and theoretical per¬ 
spective. Water-hammer experiments are a prototype model for many situations in industrial 
and military applications (e.g., trans-ocean pipelines and communication networks) where we 
have fluid-structure intera ctio n and a consequent propagation of shock-waves. After the pi¬ 
oneering work of KorteweglU (1878) and Joukowsky^ (1900), who modeled water-hammer 
waves by neglecting inertia and bending stiffness of the pipe, a more comprehensive inves¬ 
tigation, developed by Skalak m in the Fifties, considered inertial effects both in the pipe 
and the fluid, including longitudinal and bending stresses of the pipe. Skalak combined the 
Shell Theory for the tube deformation and an acoustic model of the fluid motion. He shows 
there is a coexistence of two waves traveling at different speeds: the precursor wave (of small 
amplitude and of speed close the sound speed of the pipe wall) and the primary wave (of 
larger amplitude and lower speed). Additionally, a simplified four-equation one-dimensional 
model is derived based on the assumption that pressure and axial velocity of the fluid are 
constant across cross-sections [14] . Later studies of Tijsseling [15]- m have regarded model¬ 
ing of isotropic thin pipes including an analysis of the effect of thickness on isotropic pipes 
based on the four-equation model HZ]. While all these papers consider the case of elastically 
isotropic pipes, the investigation of anisotropy in water-filled pipes of composite materials 
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was first obtained in [50] where stress wave propagation is investigated for a system composed 
of water-filled thin pipe with symmetric winding angles ±0. In the same geometry, a plat¬ 
form of numerical computations, based on the finite element method, was developed in [13] 
to describe the fluid-structure interaction during shock-wave loading of a water-filled carbon- 
reinforced plastic (CFRP) tube coupled with a solid-shell and a fluid solver. More complex 
situations involve systems of pipes mounted coaxially where the annular regions between the 
pipes can be filled with fluid. In this scenario, Biirmann has considered the modeling of non¬ 
stationary flow of compressible fluids in pipelines with several flow sections [5| . His approach 
consists of reducing the system of partial differential equations governing the fluid-structure 
interaction in coaxial pipes into a 1-dimensional problem by the Method of Characteristics. 
Later works have appeared on the modeling of sound dispersion in a cylindrical viscous layer 
bounded by two elastic thin-walled shells m and of the wave propagation in coaxial pipes 
filled with either fluid or a viscoelastic solid |H| . 

Motivated by the recent experimental effort of J. Shepherd’s group on the investigation 
of the water-hammer in annular geometries m, 0, m, m, we extend the modeling work 
of [17] and |20j to investigate the propagation of stress waves inside an annular geometry 
delimited by two water-filled coaxial pipes, in elastically isotropic and CFRP pipes. A pro¬ 
jectile impact causes propagation of a water pressure wave causing the deformation of the 
pipes. Positive extension in the radial direction of the outer pipe, accompanied by negative 
extension (contraction) in the radial direction of the internal pipe, causes an increase in the 
annular area thus activating the fluid-structure interaction mechanism. 

The architecture of the paper is as follows. After reviewing the work of You and Inaba 
on the modeling of elastically anisotropic pipes, we present the six-equation one-dimensional 
model (Paragraph |2. 2 [ ) that rules the fluid-solid interaction in a two-pipe system. In Section]^ 
we compare our theoretical findings with experimental data obtained during a series of water- 
hammer experiments. Finally, in the case of fiber reinforced pipes, the wave propagation and 
the computation of hoop and axial strain are described in full detail in Paragraph |3.3[ 




Figure 1: Schematic representation of coaxial thin pipes. LEFT: cross section. RIGHT: 
lateral view. Notice that here we are referring to the case of CFRP pipes with winding 
angles ±0. 


2 Thin pipes modeling 

2.1 One-dimensional fluid-strnctnre modeling 

According to the technique of Tijsseling one-dimensional governing equations for the 
liquid and the pipes can be obtained upon averaging out the standard balance laws in the 
radial direction. By adopting a cylindrical coordinates system, this approach is based upon 
the assumption that the behavior of water velocity and pressure depend only on the spatial 
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Table 1: Notation. Here and in what follows subscript i is either set to be equal to 1 (in the 
case in which we refer to the internal pipe) or 2 (external pipe). 


C, 

stiffness matrix (x,y and z coor¬ 
dinates) 

^i,z 

two-dimensional displacement com¬ 
ponents of pipe i along r and z axis 

ss. 

compliance matrix (x,y and z 
coordinates) 

'Ui^z 

two-dimensional velocity compo¬ 
nents of pipe i along r and z axis 

r, ip,z 

cylindrical coordinates 


one-dimensional velocity compo¬ 
nents of pipe i along r and z axis 

t 

time 

'Ui,zO 

magnitude of iii^z 

Vf 

volume fraction (fiber) 

9 

fiber winding angle 

p{r,z,t) 

two-dimensional fluid pressure 

V 

one-dimensional axial fluid velocity 

P{z,t) 

one-dimensional fluid pressure 

Vo 

magnitude of V 

Po 

magnitude of P 

Vr,Vz 

two-dimensional fluid velocity com¬ 
ponents along r and z directions 

Pout 

Pin 

pressure outside pipe 2 
pressure inside pipe 1 

£i,li 

normal and shear strain in pipe i 

Rt 

inner radius of pipe i 

c 

wavespeed 

Cj 

thickness of pipe i 


normal and shear stress 

Pi-it 

density (pipe i) 

^i,r 5 ^i,z 

two-dimensional stress components 
along r, (p and z axis 

Pw 

density of fluid 


one-dimensional axial stress 

K 

bulk modulus of fluid 

O'i,z0 

magnitude of Wi^z 

E(1),E(3) 

effective Young’s modulus along 
transverse and longitudinal di¬ 
rections in a single ply 

Em, Ef 

Young’s modulus for matrix and 
fiber 

G 31 

effective shear modulus in a sin¬ 
gle ply 

Gm,Gf 

shear modulus for matrix and fiber 


Poisson’s ratio for matrix and 
fiber 

Pm, Pf 

density of matrix and fiber 


variable z. In what follows we define one-dimensional cross-averaged quantities and obtain 
the corresponding field equations. 


2.1.1 Governing equations for the fluid 

The balance laws in the coordinate system (r, z) for the fluid read (TS] 


dvz dj) 

(2-d) Axial motion equation: Pw—;^- + 7 — = 0, 

ot oz 

dvf dp 

(2-d) Radial motion equation: Pw^z- + vw = 0, 

at or 


(2-d) Continuity equation: 


I dp dv^ I d{rvr) _ 
K dt ^ dz ^ r dr 


Here Vz{r,z,t) and Vr(r,z,t) are, respectively, the axial and radial velocity of the fluid and 
p{r, z, t) is the pressure; K is the bulk modulus of the fluid and pw is the density of the 
fluid. We now introduce the cross-sectional averaged (one-dimensional) velocity and pressure. 
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defined respectively as 


V{z,t) := 


rR2 


tt{RI - {Ri + ei)2) 


27rr Vz{r, z, t)dr, 


P{z,t) := 


rR2 


'^{^2 ~ (^1 + 6l)^) JRi+ei 


2Trrp{r, z, t)dr. 


( 2 . 1 ) 

( 2 . 2 ) 


We are in a position to introduce the one-dimensional equations of balance for the fluid, 
which are 


(1-d) Axial motion equation 

(1-d) Radial motion equation 

dVr 


dV dP „ 


(2.3) 




Rip 


r=i ?2 


-(i?i -I- ei)^p 


r=i?i+ei 


r=R2 


(Rl-iRi+eir) 


-P = 0 


(1-d) Continuity equation 


1 dP dV 


K dt ^ dz ^ {Rl - {Ri + ei)2) 


i ?2 Vr 


— (i?i -I- Cl) ' 


lr=fii+ei 


(2.4) 


= 0. (2.5) 


We remark that Eq. ( |2.4[ ) has been obtained by multiplying the two-dimensional radial 
motion equation by 27rr^, integrating in r from Ri + ci to i ?2 and dividing by 27r(i?2 ~ {Ri + 
ei)^). Here i?i is the internal radius and ei the thickness of the internal pipe while i ?2 is the 
internal radius and 62 the thickness of the external pipe (see Figj^LEFT). Moreover, in Eq. 
(2.4) it is assumed that 


dvr „ dvr 

^~di ~ ^~dt 


= («.+e,)^ 

r=i ?2 dt 


r—Ri-\-ei 


( 2 . 6 ) 


This is consistent with the (2-d) Continuity equation under the hypothesis that K is large 
and that the axial inflow is concentrated in the central axis in the limit Ri —>■ 0 , e^/Ri —^ 0 

m 


2 . 1.2 Governing equations for the pipes 


Letting i = 1,2, the equations of Axial motion and Radial motion in the pipes in the space 
(r, z) are 


(2-d) Axial motion equation: 


Piyt 




dcri,z 

dz 


= 0 , 


(2-d) Radial motion equation: 


_ 1 d{rai^r) 

Pi,t Q, 

at r or r 


Here pi^t is the density of the pipe, z, t), ai^zi'r, z, t) and z, t) are the radial, axial 

and hoop stress respectively and Ui,z are the radial and axial velocity respectively. By 
applying the cross-sectional average technique we obtain 
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(1-d) Axial motion equation 


Pi it 




dai,z 

dz 


= 0 , 


(1-d) Radial motion equation 


where 




Pi,t 


(^Ri G-i'jCTi r RiCTi 

Ri+ei Ri 


1 


dt ei{Ri-\-Ci12) ei{Ri + eil2) (i?j-|-6^/2) 




:= 


('Ri-\-e.i 


'K{{Ri + eiY - Ri) JRi 


27rr Ui^z{r, z, t)dr, 


(2.7) 


( 2 . 8 ) 


(2.9) 


nRi-\-ei 

u,^riz,t) := ——— --- 2 Trru^^r{r,z,t)dr, 

7r((i?, + ei)2 _ Rf) 


( 2 . 10 ) 


2 rRi-\-&i 

CT^,^(z,<) := -T——- 2Trrai,z{r,z,t)dr, (2.11) 

TT[[Ri + ay - Rf) Jr. 

are respectively the one-dimensional axial velocity, radial velocity and axial stress and 

pRi-\-ei 

:= — / aiJr,z,t)dr. (2-12) 

e* Jr^ 


2.1.3 Elastic properties of pipes 

By introducing the stiffness matrix and the compliance matrix SSi := the stress- 

strain relation under the plane stress assumption reads, respectively, HOI Eqs. (13, 14)] 


( ^i^x \ / Ci^ii Ci^l3 0 \ / €i^x \ 

O'z.z 1 = I C'i,13 Ci,33 0 j I Si^z j , (213) 

Ti,zx J \ 0 0 j y '^i^zx j 

( ^i,x \ / 5^1,13 0 \ / (Ti x \ 

Si,z j = I *S'i,13 S'i 33 0 j j (Ti,2 j . 

'yi,zx / \ 0 0 Si,55 / \ / 


In the case of elastically homogenous and isotropic pipes, tensor Ci reads 


/ C^,n 

Ci,13 

0 ^ 

\ / 


izi,tEi,t/ii - ufy 


C^^3 

C'i, 33 

0 

- 

v^,tE,,t/{l - izfy 

E^,tl{l-izl) 

0 

V 0 

0 

C'i,55 y 

/ \ 

V 0 

0 

G^,t J 


(2.14) 


and, in turn 


( S^,ll 

Si, 13 

0 ) 

\ I 


^i,t/ Ei,t 

0 \ 

S^,13 

Si, 33 

0 

= 

Ri,t 


0 

V 0 

0 

‘S'i,55 j 

' \ 

0 

0 

^/G^,t J 


(2.15) 
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Figure 2: Schematic rep¬ 
resentation of CFRP lay¬ 
up structures as a com¬ 
bination of single uniaxial 
plies. Axes 1, 2 and 3 are 
the principal axes of a sin¬ 
gle ply. 


where Ei^t, and Gi^t = Ei^t/+ are the Young’s modulus, Poisson’s ratio and shear 
modulus of the material from which pipe i is made. 

For anisotropic composite (fiber-reinforced) pipes the stiffness elements Ci^ki are neces¬ 
sarily a function of the geometric and elastic properties of hbers and of matrix, including the 
fiber winding angle 9. The difficulty in describing the elastic properties of hber-reinforced 
plastic thin pipes has been studied in | 20 j under the assumption that pipes are obtained by 
rolling up a woven layer with symmetric angles ±0. Each of these layers can be considered 
as a lay-up structure of multiple plies of same thickness as shown in Figure To keep the 
notation simple, in what follows we drop the subscript i. Elastic moduli are computed as a 
function of fiber volume fraction Vf and fiber and matrix elastic coefficients 0 


E(3) = EfVf + E^{1 -Vfl ^ ^ + 

A(i) Ef 

i/gi = v^Vf + Vmil - Vf), 


(1 - Vf) 1 Vf I (1 - Vf) 

Em ’ G 31 Gf ^ Gm 

Pi,t = PfVf+Pm{l - Vf) 


(2.16) 

(2.17) 


where Em,Gm,^m and pm are, respectively, the Young’s modulus, shear modulus, Poisson’s 

ratio and density of the matrix. Then Ef,Gf,iyf and pf are, respectively, the fiber Young’s 

modulus, shear modulus, Poisson’s ratio and density of the fiber. The subscripts 3 and 1 

indicate the longitudinal and transverse direction of a single ply (see Fig. [^. The Poisson’s 

ratio V 31 is defined as the ratio of the contracted normal strain in the direction 1 to the normal 

strain in the direction 3, when a normal load is applied in the longitudinal direction. The 

stiffness matrix for the composite is given by the volumetric average of the elastic stiffness 

_ 

matrices from each single ±0 ply, denoted in what follows with C . Precisely, if all plies 
have the same thickness, the stiffness matrix for a woven layer of pairs of ±0 plies is given 
by 


C=i(C+%C-^). (2.18) 

The components of C read [20l Eq. (16)] 

= Gi/ = [cos'‘( 0 )E(i)/F;ii -I-sin'‘(0)£’(3)/F;tj] +2 8 \tV{ 0 )cos^{ 9 )[v^iE^^i)/E^ -I-2G31], 
G33 = G33 = [cos^(0)£'(3)/£'jj -I-sin^( 0 )E(i)/Ejj] -|-2 sin^( 0 ) cos^( 0 ) [t'3iE(i)/i?jj -|-2G31], 

G^f = C~! = [A(3)/At) + F;(i)/F;s - 4 G 31 - 2i/3iE(i)/£;j] sin2(0) cos2(0) + 

G 55 = G 55 = [E(i)/Eu — 2 r' 3 iE(i)/E(j-I-E( 3 )/Eu — 4 G 31 ] cos^(0) sin^(0)-I-G 31 , 
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with = 1 — Finally, the compliance elements read S 33 = 

Cii/C^, Si 3 = —C i^/Cj and Cj = C 11 C 33 — where Cki are the elements of the matrix 
C defined in (2.18). 


2.2 Six-equation model 

The axial and hoop strains in pipe i can be written as [20] 

respectively. By using the strain-displacements relations 

dui. 


^i.z — 


dz ’ 


(2.19) 


( 2 . 20 ) 


by differentiating in time and by taking the cross-sectional average, Eq. (2.19 )-LEFT becomes 

da. 


dUl,z -n 

—7\ - — *i,13- 

OZ 


dt 


tp -Q dui^-^ 
-H ^i,33 r,, 

at 


where 


a,^^{z,t) := 


27T 


f'Ri -\-&i 


+ eiY - Rf) Jui 


rai^^{r,z,t)dr 


( 2 . 21 ) 


( 2 . 22 ) 


is the one-dimensional (cross-averaged) hoop stress. 

The radial displacement equation is obtained by plugging another strain-displacement 
relation, which is, 

^ (2.23) 


into Eq. (2.19)-RIGHT yielding 


(2.24) 


The equations of fluid and pipes are coupled by boundary conditions along the interfaces. 
Indeed, at each fluid-solid interface, we equate the radial velocity and radial stress of the 
fluid with those of the solid. 


(2.25) 


As in we assume that the external and internal pressures in each pipe induce a hoop 
stress which is constant in (p. Accordingly, we have m, [H] 



II 

1 

II 

to 

^'^Ar=R2 

- D 5 

' \r—R2 ’ 

1 r—Ri-\-ei 

Ar=R-i+ei ’ 

lr=i?i+ei 


^‘^^f‘\r^R2-\-e2 

= —= const., 

'^'^Ar=R2+e2 

_ yout _ Om/s) 

II 

b 

= — P*” = const., 

'^^Ar=Ri 

= I/l" = const. {= Om/s) 




1 i??(i?i+ei)2(P-P-) RlP^^-(R,+e3)^P 


2(i?i -I- ei/2)ei 


-b 


2(i?i + ei/2)ei 


a2,ip — -- 


1 RUR 2 + e 2 )HPout - P) , RlP-{R 2 + 62 ?Po^ 


(2.26) 


2 (i?2 + 62 / 2)62 


2 (i ?2 + 62/2)62 
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By plugging Eqs. (2.25 Lines 1 and 2, RIGHT) into Eq. (2.5) we obtain 
1 dP ^ dV ^ 2 

kIH {Rl - (i?i +ei)2) L 


R 2 U2- 


r=R2 


— {Rl + ei) 


r=(Ri+ei) 


= 0. (2.27) 


Now, assuming that radial inertial forces are ignored in both fluid and pipes and that the pipes 

cross-sections remain plane for axial stretches (thus implying the independency of 

on r, especially in thin pipes) we obtain a simplified model. Upon substitution of Eqs. (2.24) 


and (2.26) into Eq. (2.27) and upon substitution of Eq. (2.26) into (2.21) (by replacing 
'^'^,Ar-Ri+ei ~ '^'^’Ar-R 2 ~ ^ 2 , 2 ) we obtain new equations 

dP dV 
dz 


TO 24 

msi 

rriQi 


W 21 


dt 


50-2.2 dai^z „ 

- 17123 = 0. 


dt 


dt 


(2.28) 


and 


where 


diii.z dP dai^z 


du 2 ,z dP da 2 ,z 

^=m„—+ S„3 — 


(2.29) 


m2r := ^ 

1 

K 


2{Ri 


Q (.^2 + 62)^ -b i?2 

*J2,11 


2 (i ?2 + 62/2)62 


I I'n I . ',2 r-o -^1 + (-^1 + 

+ r-" 2(fl,+e./2K 


17123 ■= -S'!. 
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:= S' 


i?2_(i?i+ei)2r"’“ 
2{Ri + 61 )^ 

Rl - {Rl + ei)2 

2Rl 


{Rl-{Ri+eir) 

r-c -^2 , -c (-^1 + 61)^ 
*52.1 1 -r <51,11- 

62 61 


(2.30) 


Rl - {Rl + ei)^ 

<S'l,13^1 


<S'2,i3^; 


with 

Hi 


:= -In 


Ho ■■= In 


(^+ 17 ) 


(i?i + ei)2 


2i?2 


{Rl 


■61 


1 7 L2(i?i -|- ei/2)eiJ ei(2i?i -I- ei) 2(i?i -|- ei/2)ei 
(i?2 + 62)^ - 2i?2 ^ R2 

62 ' 


i ?2 7 L 2 (i ?2 + 62 / 2)62 j 62(2772 + 62 ) 2 ( 7?2 -b 62 / 2)62 


(^) 


(2.31) 


The simplified expressions in Eqs. (2.30) and (2.31) are obtained under the assumption 


(ei/(7?i -b 61 )) 1, 62/772 ^ 1- Note that it is not possible for the terms in Eq. ((2.30) to 

become singular, since the denominator becomes zero only when the annulus of water has 
zero thickness. Summarizing, the six-equations model with the six unknowns(P, V,ai^z,Ui^z) 
for the two-pipe system read 

Fluid (axial motion - continuity equation) 


514 dP „ 


dP dV 


77124- 


dcr2,. 

dt 


- 17123- 


d(Ji^, 

dt 


= 0 , 


(2.32) 

(2.33) 


8 

















































Pipes (axial motion - 2 Lxial strain equation) 


Piyt 


dai^; 
~dt 


diii 


dz 

a n c 

— -*1.33' 


= 0 , 


= 0 . 


r\ o, *-"4,00 r\, 

az at at 

We seek solutions of ( 2.32|2.35 ) in the form of wave functions, 

P = Pof{z-ct), V = Vof{z-ct), ai^z = cfi,zof{z - ct), ii, 


(2.34) 

(2.35) 

.2 = Ui,zQf{z - ct), (2.36) 


where Po,Vo,(Ti^zO and iii^zo are magnitudes and c is the wave speed. Substitution of (2.36) 
into (2.32|[235) leads to six linear homogeneous equations which we write in a compact form 


77711 

—C 

0 

0 

0 

0 


f Po \ 


f 0 \ 

-C77721 

1 

C77723 

— C77724 

0 

0 


Vo 


0 

0 

0 

-mss 

0 

—c 

0 


Cl,z0 


0 

0 

0 

0 

— 77744 

0 

—c 


C2,z0 


0 

C77751 

0 

C77753 

0 

1 

0 


Ul,z0 


0 

C77761 

0 

0 

C77764 

0 

1 


\ U2,z0 ) 


1 0 / 


(2.37) 


M 


To keep a uniform notation we have defined mu = 77133 = Pi t j ^ 7.44 = P 2 ], 77153 = S' 1^33 

and 77764 = S' 2 , 33 . Existence of non-trivial solutions to |2.37 ) requires the determinant of M 
to be zero, yielding, in turn, the following dispersion relation: 


(777247775377761 - 777237775177764 “ 777217775377764) C® + 

(777217773377764 + 777217774477753 + 777237774477751 - 777247773377761 + 777117775377764)c"‘- 
(-777217773377744 - 7771177733 77764 “ 777ii7774477753) + 777 117773377744 = 0. 


(2.38) 


Natur al frequencies of the system Ck, with fc = 1,..., 6 , are the roots of (2.38). In general, Eq. 


(2.38) has to be solved by means of numerical methods. However, if pi^t = P 2 ,t, iSius = *S' 2 ,i 3 
and S'i ^33 = S' 2,33 (e.g., if the pipes are composed of a matrix and hber with the same 


volume fraction and elastic properties) we can find exact solutions of Eq. (2.38) analytically. 
Indeed, if we define p := 77733/77753 = (/ 0 i,t<S'i, 33 )“\ g := m ii /77721 = (pu, 7772 i)“^ and <5 := 
(7772477761 — 7772377751 ) /(7775377721 ), with (5 > 0 then Eq. (2.38) I'eads 


- (c^ - -q)+ -p)=0, 


(2.39) 


with roots 


Cl = Vp, 


C2 = 


C3 = 


(p + q) 1 

1 {p + qY 

4pg(l - (5) 

^2(1-5) ) 

/ 4(1-(5)2 

4(1-(5)2 

{p + q) 

/ {p + qY 

4pg(l - (5) 

\2{l-S) \ 

/ 4(1-(5)2 

4(1-(5)2 


(2.40) 


C 4 = —Cl, C 5 = —C 2 , C 6 = —C 3 = —Cw Here ci,C 2 and Cw are positive (forward traveling) 
while C 4 , C 5 and C 6 are negative (backward traveling) wave speeds. Since c^, is smaller than 
Cl and C 2 , we refer to it as the speed of the primary wave. Accordingly, we call ci and C 2 the 
speeds of the precursor waves related to pipe 1 and 2 , respectively. 
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2.2.1 Reconstruction of the physical quantities 


We are now able to recover the m echan ical strain in the hoop and axial directions as functions 
of P, V, ai^z, Ui^z- Thanks to Eq. (2.19 l-RIGHT we can write the one-dimensional hoop strain 
as follows 


1 




^i^hoop — — 


2'!T{Ri -(- eil2)ei 


2Trrei^^dr = -I- Si^iscr^ 


where ai^^p has been obtained in (2.26) (with Pin = Pout = 0)- In turn, we have 
^i^hoop — Si^iiHiPof (^Z ct) p ct). 


(2.41) 


The cross-sectional averaged axial strain can be obtained by Eqs. (2.20) and (2.35) 

pRi-\-Ci 


1 


^i.ax — — 


^TvrSi zdr = 


dUn: 


Ui zO 


f{z - ct) = 


2'K{Ri + ei/2)ei Jji. ^ 

Si,i 3 HiPof{z - ct) + Si^33(Ji^zof{z - ct). (2.42) 

Notice that here ut^z = i—Ui,zo/c)F{z — ct) follows from integrating the last equation in 
(2.36) over time with F' = f. Then, from the system of equations (2.37) we can easily derive 
the following relations 


'^i,z0 — 




^i,zO — ^Pi,t 


(~1 + C^ Si^33Pi^t) 
and, by taking c= c^ = Cw, 

Po — Cu: Pw bo. 


cSi^l^FliPo 

(~1 + C^ Si^ssPi^t) 


(2.43) 


(2.44) 


Thanks to Eq. (2.44) we can express averaged hoop and axial stress dependent on either the 


fluid velocity Vq or, upon inversion of Eq. (|2.44|), on the fluid pressure Pq. By plugging Eqs. 
md ( 

^i.hoop — P * 5^,13 CzoPi^t 


(2.43)-RIGHT and (2.44) into Eq. (2.41) and (2.42) with c = Cw vie obtain 

CwFi^l3 Fti 


^i.ax — 


(-1 -k C^5'i,33/9i,t) 

SipsFli 


j ^Vq/*(Z Cu]t)^ 


^i,zO — 


( 1 “1“ ) 
2 


( 1 “1“ 


^wPw^ofi,^ 

Cw Pw ^0 ’ 


(2.45) 

(2.46) 


3 Water-hammer experiments 

Armed with the set of analytic expressions from Section we turn now to the simulation 
of experimental measurements for the water-hammer experiment for a set of pipes including 
homogeneous (isotropic) metal pipes and fiber-reinforced (anisotropic) pipes. 


3.1 Experimental setup 

The propagation of waves in the annular space between two pipes is studied experimentally 
using the apparatus shown in Fig. For all experiments, the outer pipe is a thick-walled 
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Pressure Transducer Ports 



Support Plate 


Figure 3: Diagram of apparatus used to measure wave propagation speeds (rotated 90° 
counterclockwise). “Electrical feedthrough” provides a sealed connection between the leads 
of strain gauges (attached to specimen tube) and the data acquisition electronics external to 
the setup. The web version of this article contains the above plot figure in color. 


cylindrical vessel made from 4140 high strength steel with an inner radius i ?2 = 38.1 mm, 
wall thickness 62 = 25.4 mm, and length 0.97 m. Tubes of various sizes can be mounted 
concentrically inside of this vessel; these tubes are held in place at both ends by polycarbonate 
plugs and sealed with gland seals. The plug at the right end of Fig. [^is fixed to the base 
of the setup while the plug at the left end is fixed to a “support plate” (Fig. inset) which 
features four holes that allow pressure waves to pass freely while still providing support for 
the inner tube. 

The annular space between the two pipes is filled with distilled water while the sealed 
cavity inside of the inner pipe remains filled with air at ambient pressure. The water inside 
the annular cavity is also at atmospheric pressure at the beginning of the test. An aluminum 
“buffer” is then inserted into the top of the apparatus, and any air that remains in the system 
is vented through a hole in the buffer that is later sealed. The buffer is composed of a 125 mm 
long aluminum cylinder capped by a 25 mm thick steel striker plate which prevents damage 
of the buffer during projectile impact; this striker plate is bolted to the aluminum buffer 
using eight 1/4-20 machine screws. 

Next a projectile is fired from a gas gun into the buffer, and the stress wave that develops in 
the buffer is transmitted into the water as a shock wave which travels along the annular space 
between the inner and outer pipes. Examples of pressure traces plotted on an a; — t diagram 
are shown in Fig. they demonstrate a sharp shock wave followed by an approximately 
exponential decay. Slight steepening of the wavefront is visible as the wave progresses. The 
high frequency fluctuations in the pressure signal are the result of axisymmetric vibrations 
of the inner pipe, as was confirmed by hoop strain measurements (not shown) at multiple 
locations around the circumference. 

The axial propagation speed of the pressure wave is measured using a row of six pressure 
transducers, PCB model 113A23, which have a response time less than 1 ^s, resonant fre¬ 
quency above 500 kHz, and are sampled at 1 MHz. These transducers are mounted flush with 
the inner surface of the outer pipe. The response of the tube is recorded using bonded strain 
gauges oriented in the hoop direction, which are coated with a compliant sealant (Vishay 
PG, M-Coat D) to avoid electrical interference from the water. Signals are amplified using 
Vishay 2310B signal conditioners and digitized at 1 MHz. 

The speed of the pressure wave is measured by marking the time of arrival of the pressure 
wave at each transducer and fitting a line to the data as shown in Fig.|^ the slope of this line 
is the wave speed. The time of arrival of the pressure wave at each transducer is determined 
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Figure 4: Example of experimental pressure traces used to determine wave propagation speed. 


as the first instant at which the pressure exceeds a chosen threshold. Because the shock wave 
steepens as it progresses, the wave speed depends slightly on chosen threshold level. For the 
results reported in this paper, wave speeds were calculated using both 50% and 80% of the 
maximum pressure of the incident wave as threshold values, and in every case the results 
differed by less than 5%. The lower threshold value of 50% was chosen because it is larger 
than the amplitude of the precursor wave, which ensures that the measured primary wave 
speed is not influenced by that of the precursor wave. 

For each specimen tube, between 3 and 8 shots were conducted. After each shot, wave 
speeds were calculated using both 50% and 80% of the maximum pressure as the threshold 
described above, and the average of these two values was taken as the measured wave speed 
for the shot. Finally, the average wave speed over all 3-8 shots was calculated, and this 
number was taken to be the wave speed associated with the specimen. The uncertainty in 
the measured wave speed was calculated as twice the standard deviation of the measured 
wave speeds over all 3-8 shots; this uncertainty is typically less than 10%. The difference in 
measured wave speed between shots depends mainly on the quality of the impact between the 
projectile and buffer. Slight non-normal impact is unavoidable and leads to slight differences 
in the measured velocity. 

Further details about the experimental setup are recorded in Refs. 011]. It may be noted 
that the experiment was designed for the study of plastic deformation and buckling of tubes 
loaded by dynamic external pressure, which is the reason for the rather high pressure of the 
shock wave (on the order of 1-10 MPa). For all results reported in this paper, the pressure 
and total impulse of the pressure load were low enough to prevent plastic deformation of both 
the inner and outer pipes, as was verified by hoop strain measurements. Although the hoop 
strain measurements indicated slight non-axisymmetric deformation of the inner pipe (elastic 
buckling) in most shots, these effects do not appear to significantly influence the propagation 
of pressure waves so long as the buckles remain elastic (see Ref. 0). 
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Table 2: Parameters of the large 4140 steel tube (referred as pipe 2 in Sectionand water. 


4140 Steel (pipe 2) 

Young’s modulus 

E2,t 

198 GPa 

Inner radius 

i?2 

38.1 mm 

Density 

P2,t 

7800 kg/m3 

Thickness 

62 

25.4 mm 

Poisson’s ratio 

V2,t 

0.3 




Water 

Bulk modulus 

K 

2.14 GPa 




Density 

Pw 

999 kg/m3 





Table 3: Parameters of stainless steel and aluminum. 


Aluminum (pipe 1) 



Stainless steel (pipe 1) 



Young’s modulus 

Ei^t 

68.9 GPa 

Young’s modulus 

Ei,t 

198 GPa 

Density 

Pi,t 

2700 kg/m^ 

Density 

Pi,t 

8040 kg/m^ 

Poisson’s ratio 


0.33 

Poisson’s ratio 

Vl,t 

0.29 


3.2 Results and discussion 


We now address our experimental measurements of the speed of the primary wave, the fluid 
pressure and the hoop strain in the internal pipe against the prediction based on our modeling 
work. The geometrical and physical properties of the external pipe are reported in Table 
The smaller interior tubes are mounted concentrically and are made from either aluminum, 
stainless steel (with properties listed in Tables and [4|) o r a carbon fiber-epoxy resin matrix 
composite. The elastic coefficients matrices and SSi for both the internal (i = 1) and 
external {i = 2) pipe are computed as in Paragraph |2.1.3| with the only exception being the 
CFRP pipe for which a clarification of the method is given in the Remark below. Notice 
that since the physical properties of the external and internal pipes are different, and in 
particular Ci ^ C 2 , computation of the primary wave speed Cw needs to be accomplished 
by solving Eq. (2.38) numerically. Results displayed in Table show a good match between 
the experimental data and the predicted values. However, there are some slight differences; 
in particular, the measured wave speeds are consistently greater than those predicted by the 
model. One possible explanation is that the steepening wavefront of the pressure wave biases 
the experimental measurements toward higher wave speeds. 

Comparisons of experimental measurements and predictions for hoop and axial strain in 
the internal pipe, given in Table show that the calculated data match with experimental 
results reasonably well. In this table the measured data points are the peak values of the 
pressure and hoop strain, which were determined after applying a 50 kHz low-pass filter to the 
recorded signals. The computed values of these physical variables are based on Cw since the 
primary wave generates larger magnitudes in both axial and hoop strain than the precursor 
waves. Differences between the computed and measured values are the result of modeling 
idealizations that are not completely satisfied in the experiment. For instance, radial inertia 
of the fluid and pipes, strain rate effects, axial bending of the pipe in the vicinity of the 
wavefront, non-planar features of the pressure wavefront, and transient phenomena, may 
well contribute to the differences between measurement and computation. 

Although only direct measurements for the fluid pressure are available, we are able to 
report estimated data for the fluid velocity as well. Since a direct measurement of the fluid 
velocity is not available, the fluid velocity behind the wavefront is estimated as being equal 
to the initial velocity of the buffer (see Fig. 3). The initial buffer velocity was estimated 
by equating the momentum of the projectile prior to impact with that of the buffer after 
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impact. In other experiments using the same experimental setup, high speed video of the 
projectile-buffer impact has shown that this estimate of the buffer velocity is typically within 
10 -20% of its actual velocity. 


Table 4: Geometrical parameters of pipe 1 (six different cases). f = aluminum, | = stainless 
steel, * = carbon-epoxy composite. 


Tube ID 

Ri [m] 

Cl [mm] 

Speed of primary 

wave Cuj [m/s] 

Measurement 

Computation 

2 ff 

0.0199 

1.47 

1245 ± 14.72 

1207 

34t 

0.0218 

0.89 

1078 ± 123.50 

1054 

351- 

0.0154 

0.89 

1309 ± 36.50 

1281 

381' 

0.0154 

0.89 

1312 ±46.76 

1281 

3fp 

0.0154 

0.89 

1379 ± 95.92 

1369 

40* 

0.019 

1.45 

1157 ±38.03 

1100 


Table 5: Comparison of measurements and predictions of maximum hoop strain accompa¬ 
nying the water-hammer wave (pipe 1) and fluid velocity. Here the fluid velocity Vq was not 


measured, but instead was estimated at the post-processing level. Variables with a super¬ 
script ^ have been computed using the water pressure Pq determined from the experiments 
and by using Eqs. (12.441) and (12. 451). 


Experiments 


Computations 


Shoot ID 

Po 

[MPa] 

^l^hoop 

[mstr] 

[m/s] 

^l,hoop 

[mstr] 

[m/s] 

100 (tube ID 27) 

3.10 

-0.61 

2.86 

-0.6581 

2.57 

106 (tube ID 27) 

6.44 

-1.23 

5.69 

-1.3672 

5.34 

137 (tube ID 34) 

2.64 

-0.93 

2.83 

-0.9816 

2.51 

151 (tube ID 35) 

4.98 

-1.37 

3.97 

-1.3325 

3.89 

200 (tube ID 38) 

3.05 

-0.55 

3.14 

-0.8161 

2.38 

159 (tube ID 36) 

7.92 

-0.61 

5.22 

-0.7054 

5.76 

160 (tube ID 36) 

11.10 

-0.76 

6.19 

-0.9886 

8.07 

398 (tube ID 40) 

3.12 

-1.56 

- 

-1.3 

- 


Remark 1. While our analysis is valid for pipes of various thicknesses and elastic moduli, 
experimental results are available only for a system where thin internal pipes are coupled to 
a relatively thick and stiff outer pipe. In this regime, the external pipe is almost rigid and the 
main contributions to the fluid-structure interaction are due to the internal pipe. To clarify 
this point, we analyze the dependence of the primary wave speed on certain parameters of 
the external pipe. As an example, we limit our analysis to systems with an internal pipe 
made of stainless steel with parameters reported in Table 3. Similar results can be obtained 
by considering an aluminum or composite pipe. 

In Fig. I^LEFT we report the computation of the primary wave speed as a function 
of the thickness of the external pipe and considering the Young’s modulus of pipe 2 as a 
parameter. Therefore the continuous curve corresponds to 4140 steel while the other curves 
correspond to a material with increasing stiffness including the limit case of an ideally rigid 
material. Interestingly, the computed value of the primary wave speed for a 4140 steel pipe 
with internal radius i?i = 38.1 mm and thickness of ei = 25.4 mm {cw = 1369 m/s, Table|^ 
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Figure 5: Comparison of computed primary wave speed c^, and experimental measurements. 
LEFT: aluminum, tube IDs 34 and 35 (prediction reported as continuous line); stainless steel, 
tube ID 36 (prediction reported as dashed line). RIGHT: carbon-epoxy composite, tube ID 
40 (prediction reported as a continuous line); aluminum, tube ID 27 (prediction reported as 
dashed line) and ID 38 (dash-dotted line). The web version of this article contains the above 
plot figures in color. 




Figure 6 : LEFT: Computation of the primary wave speed c^, as a function of the thickness 
of the outer tube. Here various plots correspond to increasing values of the Young’s modulus 
of pipe 2 including the asymptotic case of an ideally rigid outer pipe. The experimental 
measurement is reported for tube ID 36. RIGHT: Computation of the primary wave speed 
Cw for different values of the thickness of the outer tube. The experimental result available 
for tube ID 36 is to be compared with the curve obtained for 62 = 25.4 mm. 


is reasonably close to asymptotic value obtained for a perfectly rigid external pipe (c^ = 1404 
m/s) thus confirming that the external pipe is, with a good approximation, almost rigid. 

Fig. [^RIGHT shows the dependence of on 62 to clarify the effect of a softer outer 
pipe. Again, we consider an internal pipe made of stainless steel and with ei = 0.89 mm 
while the external pipe is in 4140 steel and with fixed R 2 = 38.1 mm. The plot of the wave 
speed, as a function of Ri, generalizes the one in Fig. [^LEFT as the thickness of the outer 
pipe is now regarded as a parameter. Indeed, for 62 = 25.4 mm we obtain the dashed curve 
plotted in Fig. [^LEFT for which a comparison with an experimental result is available. For 
62 = 2.54 mm we have R^j&i = 15 which may well be regarded in the regime of thin pipes. 
As 62 decreases we obtain decreasing profiles and values of the wave speed and, by Eq. (2.44), 
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of the fluid pressure as well. Indeed, softer pipes undergo large radial deformations causing 
an increase of the annular area and consequently a drop of pressure. 


Remark 2. The computations of the elastic stiffness for carbon-epoxy composite tubes 
deserves a special comment. When we have more fibers in one direction than in another we 
must adapt the method used in Paragraph |2.1.3[ Indeed, since tube 40 is a lay-up structure of 
layers containing fibers in both axial {6 = 0°) and transverse {6 = 90°) direction, its stiffness 

matrix, denoted by Cj^ , is determined by a volumetric average of the two stiffness matrices 

_ |_^ 

obtained by plugging 9 = 0° or 6 = 90° into C . Although the exact proportion of plies with 
fibers in either direction is unknown, we have been able to estimate the elastic coefficients 

-C6 -0*^ 

by simply assuming that the known coefficient Ci 33 is given by a linear combination of 
and Cgg [TU]. Indeed, by solving the system of equations 

+ C“33 = e^33 + (l-e)C3r, (3-1) 


we are able to compute the two unknown variables C 3 n and The latter, is the relative 
amount of fibers in the axial direction, a non-dimensional parameter which takes into account 
both the percentage of plies with fibers in the axial direction and their thickness. For pipes 
composed of uniaxial layers with all fibers in direction 0 = 0 ° we have ^ = 1 , while in the 
dual case of fiber-reinforced pipes with winding angle 6 = 90° we have ^ = 0. Notice that the 
elastic coefficient for CFRP pipes in Paragraph |2.1.3| have been computed for ^ = 0.5 since we 
have the same amount of plies with winding angle -\-9 as those with winding angle —9 and all 
the plies have the same thickness. By plugging C 33 = = 142 GPa and = C 33 = 9 

GPa [24] and = 117 GPa [23] yields « 34 GPa and ^ « 0.8. Finally, the compu¬ 

tation of the density of the composite and of the Poisson’s ratio follow as in Eq. (2.17) by 
assuming pf = 1770 kg/m^, I'f = 0.2 (carbon fiber) and pm = 1208 kg/m^, ly,. 


resin) [^ with fiber volume fraction Vf = 2f 3 


= 0.39 (epoxy 
yielding r'gi = 0.26 and pi^t = 1583 kg/m^. 


Investigation of the pressure waves and mechanical strain profiles becomes particularly 
transparent when the system is composed of two coaxial carbon-reinforced pipes and when 
fibers and matrix have the same physical properties in both pipes. Even though for this case 
experimental results are not available, in the following we report the analysis for the readers 
convenience thus generalizing the one-pipe modeling work of m- 


3.3 Fiber-reinforced plastic pipes 


The study of the fluid-structure interaction in anisotropic structures, considered for the first 
time in this paper with the modeling tube ID 40 in Paragraph |3.2[ is now analyzed in full 
detail for a system of two coaxial water-filled CFRP pipes composed of pairs of ±0 plies 
of same thickness. To keep our analysis simpler, we assume that the physical properties of 
matrix and fibers (including the winding angle) a re iden tical to one another in both pipes (see 
Tablej^ yielding Ci^ki = C 2 ,fez- As in Paragraph 2.1.3 we drop the subscript i in the stiffness 
and compliance coefficients. The general case of pipes with different properties, including 
different winding angles 0 i ^ 02 , can be studied with similar techniques and is left to the 
brave reader. 

The plots of the stiffness and compliance elements Cki and Ski as a function of the winding 
angle 0 are contained in j^Dj. Here it is enough to recall that at 0 = 0°, the direction of the 
fibers coincides with the axial direction of the pipes and consequently the axial stiffness C 33 
has a maximum while the hoop stiffness Cn has a minimum. As 0 increases, Cn increases 
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while C 33 decreases due to the fiber reinforcement in the hoop direction. Conversely, as 9 
increases, S'!! decreases while S'33 increases. Then, the absolute values of the coupling terms 
C 13 and 5 'i 3 have a maximum at 6 = 45° and are minimized at 9 = 0° and 90°. 


Table 6 : Geometrical and physical data for coaxial CFRP pipes, pO]. 


Carbon fiber 



Epoxy resin (matrix) 



Young’s modulus 

Ef 

238 GPa 

Young’s modulus 

Em 

2.83 GPa 

Density 

Pf 

1770 kg/m3 

Density 

Pm 

1208 kg/m3 

Poisson’s ratio 


0.2 

Poisson’s ratio 

^m 

0.39 

Composite pipes 

Inner radius (pipe I) 

Ri 

19.15 mm 

Inner radius (pipe 2) 

R 2 

54 mm 

Thickness (pipe 1, 2) 

61=62 

1.66 mm 

Piber volume fraction 

Vf 

0.7 


The calculated primary and precursor waves are shown in Fig. [^as functions of 9. Similar to 
the one-pipe scenario of [ID] , as 9 increases the result is a reinforcement of the hoop stiffness 
and an increase in the primary wave speed which corresponds to the breathing mode of the 
pipes. Conversely, the speed of the precursor waves (corresponding to a longitudinal mode 
of the pipes) diminishes due to decreased axial stiffness. We remark that the computation 
of the physical variables is based on the incident wave c = because its magnitude is much 
larger than those of precursor waves. In this case, by Eq. (2.44) we have Pq/Vq = and 
therefore the graph of the fluid pressure is not reported. Precursor wave speeds ci and C 2 
are very similar when the parameter 6 is very small which happens for either small or large 
9. This is consistent with Eq. (2.40) where for <5 = 0 we have ci = C 2 . It is possible to 
prove that this corresponds to the cases for which the coupling stiffness S '13 as a function of 
9 is minimal. On the other hand, it follows that for intermediate values of 9 wavespeed Ci 
becomes slower while 5 becomes larger and is maximized at 0 « 43° when also S '13 is large. 
Proof of the dependence of 5 on the coupling term S '13 requires a more detailed asymptotic 
analysis and it is therefore left to a forthcoming paper. Plots of the hoop and axial strain 
in pipes I and 2 calculated from Eqs. (2.45) and (2.46) are reported in Fig. As the 
impulsive impact by the projectile generates a positive water pressure, pipe 2 undergoes a 
positive expansion in the radial direction (hoop strain is positive) while pipe 1 is contracted 
in the radial direction (hoop strain is negative). Note that the radial expansion of pipe 
2 accompanies the contraction in the axial direction while the radial contraction of pipe 1 
accompanies the expansion in the axial direction. This explains why hoop and axial strain 
have opposite signs. The hoop strain is essentially determined by the hoop compliance An 
and therefore the absolute values of the hoop strains in both pipe I and 2 are large for small 
9 and are decreasing for increasing 9. Then, at a first order of approximation, the axial strain 
is mainly determined by the coupling term A 13 uni Sect. 5] and therefore the axial strains 
in both pipe I and 2 are maximized (in absolute value) for 9 « 52°. 


4 Summary and future perspectives 

We have investigated the propagation of stress waves in water-filled pipes in an annular ge¬ 
ometry. A six-equation model that describes the fluid-structure interaction has been derived 
and adapted to both elastically isotropic and anisotropic (fiber-reinforced) pipes. The nat¬ 
ural frequencies of the system (eigenvalues) and the amplitude of the pressure and velocity 
of the fluid, along with the mechanical strains and stresses in the pipes (eigenvectors) have 
been computed and compared with experimental data of water-hammer tests. It is observed 
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Figure 7: LEFT: plot of the primary wave speed Cw, RIGHT: plot of the precursor wave 
speeds ci and C 2 as a function of the winding angle 9 and of the coefficient S (non-dimensional, 
multiplied by 2000) as a function of 9. The web version of this article contains the above 
plot figure in color. 


X 10 " X 10“^ 



6 [deg] 9 [deg] 


Figure 8: LEFT: plot of the maximum hoop strains accompanying the water-hammer wave 
normalized by Vb = 1 m/s as a function of the winding angle 9 (Eqs. (2.45)). RIGHT: plot 
of the maximum axial strains accompanying the water-hammer wave normalized by Vq = 1 
m/s as a function of the winding angle 9 (Eq. (2.46l). 


that the projectile impact causes a positive expansion of the external pipe in the radial di¬ 
rection and a contraction in the axial direction (Poisson’s effect). Vice versa, the internal 
pipe is contracted in the radial direction and expanded axially. In the last section of the 
paper, which is a benchmark for future experimental investigation, we have analyzed in full 
detail the propagation of waves in CFRP pipes with a special emphasis on the influence of 
the winding angle on the wave speeds and the axial and hoop strains in the pipes. Most 
interestingly, we found that the speed of the primary wave (breathing mode) increases with 
the increasing winding angle due to increasing hoop stiffness in both pipes. This is in agree¬ 
ment with the one-pipe model and analysis of [20] • Gonversely, the profile of the speed of the 
precursor waves (longitudinal modes) is large when the winding angle is small (and therefore 
the axial stiffness is large) and it decreases with increasing winding angle. Additionally, we 
have observed that the two precursor waves travel at almost the same speed when the fiber 
winding angle is equal to either 0° or 90° while a separation of the velocities is observed for 
9 in between these values, a phenomenon which is the object of further analysis. 
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